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Abstract 

We construct a solution to a 2 x 2 strictly hyperbolic system of 
conservation laws, showing that the Godunov scheme |12| can produce 
an arbitrarily large amount of oscillations. This happens when the 
speed of a shock is close to rational, inducing a resonance with the grid. 
Differently from the Glimm scheme or the vanishing viscosity method, 
for systems of conservation laws our counterexample indicates that no 
a priori BV bounds or L 1 stability estimates can in general be valid 
for finite difference schemes. 



1 Introduction 

Consider a strictly hyperbolic m x m system of conservation laws in one 
space dimension: 

«t + /(«)x = 0. (1.1) 

For initial data with small total variation, the existence of a unique entropy 
weak solution is well known |11| . [7j, [3]- A closely related question is the 
stability and convergence of various types of approximate solutions. For 
vanishing viscosity approximations 

u t + f(u) x = eu xx , (1.2) 
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uniform BV bounds, stability and convergence as e — > were recently es- 
tablished in [I]. Assuming that all the eigenvalues of the Jacobian matrix 
Df(u) are strictly positive, similar results are also proved in [2] for solutions 
constructed by the semidiscrete (upwind) Godunov scheme 

^«j(t) + ^[/(«i(t))-/(«j_i(t))] =0, Uj (t) =u(t, jAx). (1.3) 

Recently, Bianchini has succeeded in proving the same type of BV bounds 
for the Jin-Xin relaxation model [3]. 

In the present paper we study the case of fully discrete schemes, where 
the derivatives with respect to both time and space are replaced by finite 
differences. In sharp contrast with the previous situations, we show that 
the total variation of an approximate solution constructed by the Godunov 
scheme can become arbitrarily large. This typically happens when the so- 
lution contains a shock along a curve x = j(t), whose speed remains for a 
long time very close (but not exactly equal) to a given rational, for example 

1 Ax 
^2 AT 

In this case, a resonance occurs and a substantial amount of downstream 
oscillations is observed in the numerically computed solution. For general 
m x m systems, our counterexample indicates that a priori BV bounds and 
uniform stability estimates cannot hold for solutions generated by finite dif- 
ference schemes. It leaves open the possibility that these difference schemes 
still converge to the unique entropy weak solution. In any case, a rigorous 
proof of this convergence cannot rely on the same arguments as in , based 
on uniform BV bounds and Helly's compactness theorem. 

We recall that, for the 2x2 system of isentropic gas dynamics, the 
convergence of Lax-Friedrichs approximations was proved in [0] , within the 
framework of compensated compactness. Further results were obtained for 
straight line systems, where all the Rankine-Hugoniot curves are straight 
lines. For these systems, uniform BV bounds, stability and convergence 
of Godunov and Lax-Friedrichs approximations were established in |17j . |Hj 
and [U- The analysis relies on the fact that, due to the very particular 
geometry, the interaction of waves of the same family does not generate 
additional oscillations. 

It is interesting to understand why the arguments in [3] or |2j break 
down, when applied to fully discrete schemes. A key ingredient in the anal- 
ysis of vanishing viscosity approximations is the local decomposition of a 
viscous solution in terms of traveling waves. To achieve a good control the 
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oscillations produced by interactions of waves of a same family, it is essential 
that the center manifold of traveling profiles has a certain degree of smooth- 
ness. This is precisely what fails in the case of fully discrete schemes. As 
remarked by Serre [2U], for general hyperbolic systems the discrete shock 
profiles cannot depend continuously on the speed a, in the BV norm. In 
the related paper we constructed an explicit example showing how this 
happens. 

Our basic example is provided by a 2 x 2 system in triangular form 

u t + f(u) x = 0, (1.4) 
v t + Xv x + g(u) x = 0. (1.5) 

This system is strictly hyperbolic provided that f'(u) > A for all u. Choosing 
mesh sizes At = Ax = 1, the Godunov (upwind) scheme takes the form 

< +1 = «?"[/(«?)-/(«?-!)]. (1-6) 
= A«y_! + (1 - A)i£ - [g(u]) - g^)] . (1.7) 

We shall assume that 

< A < f'(u) < 1 for all m£R (1.8) 

so that the usual linearized stability conditions are satisfied. Thanks to the 
triangular form of the system (|1.4[) - (|1.5|) . the exact solution of a Cauchy 
problem can be computed explicitly. Indeed, one first solves the scalar 
equation for u, say by the method of characteristics. Then the function 
g(u) x is plugged as a source term into the second equation, which is linear 
in v. 

More specifically, we shall consider a solution of ()1.4l) - (|1.5j) where the 
first component u contains a single shock, located along the curve x = 7(t). 
One can arrange things so that the speed of the shock varies slowly in 
time, remaining close to some rational number, say j(t) ~ 1/2. The second 
component v will then satisfy a linear transport equation with a source 
located along a smooth curve, and a priori bounds on its total variation can 
be easily given. 

For the corresponding approximate solution generated by the Godunov 
scheme (jl.6|) - ()1.7jl . however, things turn out to be quite different. Because 
of the discretization, in the equation for the second component the source 
is not located along a smooth curve, but sampled at grid points. Since the 
speed of the shock is close to rational, this source "resonates" with the grid, 
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producing an arbitrarily large amount of downstream oscillations as time 
progresses. 

The plan of the paper is as follows. Section El contains a preliminary 
analysis of the heat equation with a moving source. We show that if the 
point sources are located at grid points and have an average speed close to 
rational, then the solution will contain downstream oscillations. By study- 
ing the strength and the location of these tail oscillations one gets a basic 
understanding of what happens for solutions generated by discrete schemes. 

In Sections El - El we carry out a detailed construction of a Godunov 
approximate solution for the system l|1.4|) - l|1.5|) . showing that the total vari- 
ation can become arbitrarily large. Choosing the flux 

f(u) = ln((l-u) + ue u ), (1.9) 

one can perform a nonlinear transformation introduced by Lax [T3|, |15j and 
explicitly compute the discrete Godunov solution -u™, in terms of binomial 
coefficients. In Section El we summarize the main features of our example, 
and discuss its significance toward a rigorous theory of discrete schemes 
for hyperbolic conservation laws. For readers' convenience, some results on 
the approximation of the binomial distribution in terms of heat kernels are 
recalled in an appendix. 



2 The heat equation with a moving source 

As a motivation for the following analysis, consider a solution of the finite 
difference scheme 



n+l 



q - [/(«?) - /(«£.!)] 



(2.1) 



in the form of a discrete shock with speed a, say 

u] = <t>(j - crn) . 

For a scalar conservation law, the existence of discrete traveling profiles was 
proved in ^3]. We assume that the flux / is strictly convex and that the 
shock profile satisfies 



lim <p(s) > lim i 

s^— oo s— >+oo 



u 



Inserting this solution as a source in the second equation (|1.7|) . we obtain 

vf l = At^ + (1 - \)v] - [g(u]) - g^)] 

= Ai£_! + (1 - \)v? - \g(<f>(j - an)) - g(cf>(j - 1 - an))] . (2.2) 
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Figure 1: Discrete sources at integer points. 



We can choose a function g such that 

g (u) = if \u 
g (u) > if \u 

for some e > small. With this choice, the source terms on the right hand 
side of (|2.2j) will vanish outside a thin strip centered around the shock. 
Indeed, if the shock is located along the line x = at, a fairly good approxi- 
mation is 

{f + >-*-> (2.3, 

Here [s] denotes the largest integer < s. To understand the oscillations in 
(J2.2|) determined by the presence of these moving sources, we first study some 
model problems related to the standard heat equation. The approximation 
of a finite difference scheme by means of a second order diffusion equation is 
indeed a standard tool of analysis (see p. 117). In this section we review 
some calculations that highlight the mechanism responsible for generating 
variation in solutions to fully discrete schemes. For a detailed treatment 
we refer to the paper [T. Rigorous estimates will then be worked out in 
Sections 01 - El 



u + +u 
2 



> e 
< e, 
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First, consider the easier case of the heat equation with point sources 
located on a discrete set of points P n = (n, an), with n integer (the white 
circles in Figure 1) 

vt - v xx = 5 n ^ n . (2.4) 

We assume that a > and consider a solution of (|2,4|) defined for t £]— oo, 0[. 
Its values at time t = are now computed as 

v(0,y) =*(y) :=^G(n, y + an) . (2.5) 

n>l 

Here 



Git, x) := — =e 



-x 2 /4t 

is the standard heat kernel. To understand how the oscillations of <&(y) 
decay as y — > — oo, we express the sum in (|2.5j) as an integral 

/•oo 

:= V G(n, y + an) = / G{t, y + at) (l + h\ (t)) dt (2.6) 

n>l ^ 

where 

^i(t) := [*I -t + 1/2. (2.7) 

By induction, we can find a sequence of periodic and uniformly bounded 
functions h m such that 

h m (t) = h m (t + 1) , J h m (t)dt = 0, —h m (t) = h m ^ 1 (t) . 

Integrating by parts we obtain 

roc / d m \ 

<%) = J G(t,y + at)(l + —h m {t))dt 

-i /•oo jm 

= _ + (-l)*»jf — G (t,y + at)h m (t)dt. 

The identities 

G(t, x) = r l ' 2 G{l, x/Vt) , G t = G xx , 

imply 

^L G(tx ) = j-fm+D/2 d m G(l, x/Vt) 

Q x m \ ' > Q x m 
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In addition we observe that, 



sup 

\ t+ y\ < lylH-l/2 



dP 



: G(t, y + at) 



0(1) ■ e Ccl 



as y —* — oo . 



for some constant c e > 0. Letting y — > — oo, for every m > 1 the above 
estimates imply 



l 
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< 



dt' 



r G(t, y + 



dt = 0(l)-y 



-m/2 



Similarly, 



■/ 



d 



771+ 1 



dt m+1 



G(t, y + at) 



dt = 0(1) -y 



-(m+l)/2 



(2.8) 



(2.9) 



Taking m = 2 in (|2,9[) one obtains the integrability of <£>', hence a bound on 
the total variation of <E>. 

Next we outline the case when the point sources are located not at 
the points P n = (n, an), but at the points with integer coordinates Q n := 
(n, [an]) (the black circles in Figure 1), 

Vt - V xx = Sn^an} ■ 

Again we consider a solution defined for t G ] — oo, 0] and a direct compu- 
tation yields 

«(0, V-l) = *(y):=Y, G ( n >V+ M • 
n>l 

Because of l|2.8jl . to determine the asymptotic behavior as y — > — oo, it 
suffices to estimate the difference 

K(y) := tf(y)-$(y) = - ^ [G(n, y + an) - G(n, y + {an}) 

n>l 

It is here that, if the speed a is close to a rational, a resonance is observed. To 
see a simple case, let a = 1 + e, with e > small. Then we can approximate 

K (y) ~ -^G x (n,y + an)(an-\an}) 

n>l 

POO 

« - / G x (t, y + at) (et - let}) dt . 
Jo 



(2.10) 
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The functions appearing in the above integration are shown in Figure 2. 
We recall that 

for every y < 0. Set y e := —e~ 2 . When y ranges within the interval 

I £ :=[y £ , y £ /2] = [- £ ~ 2 , - e~ 2 /2] , 
the integral in (|2.10j) can be of the same order of magnitude as 

/ \G x (t, y e + at)\dt > c y~ 1/2 = c e. 
Jo 

Moreover, each time that y increases by an amount Ay = the phase of 
the fractional part \ey\ — ey goes through a full cycle, hence the map 

/ G x (t, y + at) (et - let}) dt 
Jo 

oscillates by an amount > c\S. In all, we have approximately 1/2e cycles 
within the interval I e . Hence the total variation of the discrete profile \P = 
y|r(i+e) can \y e estimated as 

Tot. Var. {^( 1+£ ) ; J £ } > c 2 (2.11) 



S 




for some constant C2 > independent of e. By (|2.11j) it is clear that, as e — > 
0+, the functions \I/( 1+E ) do not form a Cauchy sequence and cannot converge 
in the space BV. The details of the preceding argument are given (for the 
Lax-Friedrichs scheme) in . Finally we consider the case where the source 
travels with a carefully chosen variable speed. In the previous computation 
the source travelled with a constant speed a = 1+e, and a significant amount 
of oscillation was observed at a distance \y e \ = e~ 2 downstream from the 
shock. To construct a solution whose oscillation becomes arbitrarily large, 
the idea is to choose a large time interval [0, T] and generate a source with 
non-constant speed, say located along a curve x = j(t) with 

1 (t) = (t-T)-2VT^t, Tr(t) = l + _J= (2.12) 

for < t < T — T 1 / 2 , (see Figure 3). Consider the profile at the terminal 
time T of a solution of 

vt + v xx = 5 n> | 7(n) ] v(0, x) = 0, (2.13) 

defined for t £ [0, T]. For a given point y the value v(T,y) is essentially 
determined by the sources located inside the parabolic region 

r„ :={(*,*) : xe[y-(T-ty +1 / 2 , y + (T-ty +1 / 2 ]}. 
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By construction, if — T/2 < y < — T 1//2 , this region contains a portion of the 
curve 7 traveling with speed 

7(t) = l + -^=»l + -i=. 

VT-t J\y\ 



This is precisely the resonant speed that can produce downstream oscilla- 
tions at a distance ~ y from the shock. As in (|2.11|) . one has 

Tot. Var. {v(T, •) ; [y, y/2}} > c 2 -T/2<y< -T 1 ' 2 . (2.14) 

Now consider the points 

y^-VT 1 ' 2 j = 0,l,...,N-l, 

where N is the largest integer < ^ log 2 T. Notice that this choice of N 
guarantees that —T/2 < yjy < . . . < y\ < y$ = — T 1 / 2 . Observe that the 
choice of the variable speed 7 at (|2.12j) is precisely what is needed in order 
to produce a uniformly positive variation on each of the intervals [yj, 2/7-1]. 
From (j2.14|) it thus follows 

N /1 \ 

Tot. Var. {v(T, •)} > ^Tot. Var. {^(r, •) ; [ Vj , y 3 /2}} > c 2 - 1 - log 2 T - 2 J 

(2.15) 

Choosing T large, we expect to find numerical solutions generated by the 
Godunov scheme whose total variation grows by an arbitrarily large amount. 
In the following sections we show that this is indeed the case. By providing 
rigorous estimates on all the approximations performed in the above formal 
analysis, we will prove that the discrete Godunov solution satisfies the same 
type of estimate as in (|2.15|) . 

3 Special solutions 
3.1 Inviscid solution 

It will be convenient to work with the particular scalar conservation law 

u t + [ln(v + »e u )] x = 0, (3.1) 

where fi G] 1/2,1 [ and v := 1 — /i. This equation will serve as the first 
equation in the system (|1.4j) - (|1.5|) . Notice that in this case the flux f(u) = 
ln(i/ + fie u ) is increasing and convex. Indeed, 

J W = — i u ' / W = 7 — ; ^2 ' 

v + jjLe u [y + \xe u y 
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so that 

/'(u)€]0,l[, /»>0. 

Note that 

/(0) = 0, /'(0) =fi> 1/2. 

Our first goal is to construct an exact solution of the inviscid equation 1)3. 1|) 
containing exactly one shock, located along a curve j(t) which travels with 
a speed close to a rational number Ao > A, where A is as in 1)1.5)1 , The role 
of this curve will be that its spatial translates are close to the level curves 
of a certain (everywhere defined) solution of the upwind scheme (|1,6|) . This 
particular solution of the scheme will then be a good approximation of the 
exact solution of 1)3.1(1 which we construct here. 

We now fix the speed A in 1)1. 5|) to be 1/2. (This choice is not essential 
but it simplifies the calculations below which involve approximating binomial 
distributions with the heat kernel.) The shock solution of 1)1.4)1 will connect 
a variable state U-(t) to a fixed right state tt+ < u-(t), and we want this 
u-shock to travel faster than the advection speed A = 1/2 of the second 
equation. We begin by fixing the right state u+ := 0, and our solution will 
be constant equal to this value to the right of 7, i.e. 

u(t,x) = for x > 7(t). 

Next, in analogy with the analysis in Section [2J we insist that the shock 
curve 7(i) should have a speed as in (j2.12)) . For notational convenience we 
will make the explicit construction for negative times. We thus fix a large, 
positive (integer) time T, and define the curve j(t) by setting 

7 (i) = \ Q t - 2V^i, forte [-2T,-VT [. (3.2) 

We require 7 to be smooth on [—IT, 0] and such that 

7 (t) = j(-VT) = Ao + forte [-VT,0], 

see Figure 3. Here Ao £]/J,l[nQ denotes a fixed, rational speed. At every 
time t £ [ — 2T, 0] the state u-(t) immediately to the left of the shock must 
then satisfy the Rankine-Hugoniot equation 

ln(i/ + /ie"-W) 

=7() ' ( } 
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To establish the existence and properties of the function t ^ u~(t), we let 
uq denote the left state corresponding to constant speed Ao, i.e. 

— Ao , 

Uq 

and then consider the equation 

u 

A straightforward argument shows that (|3.4|) has exactly one solution u = 
y?(o") G ] 0, oo [ for each a G ] /i, 1 [, and that ip depends smoothly on a. There- 
fore, for t negative and sufficiently large, from (|3.3j) and (j3.2j) we recover a 
smooth function t u—(t), where 

u-(t) = (p(\ + 



~ Ko+ (Z^) + (Z^' W 
for t < 0. We have here assumed the Taylor expansion 
^(A + 5) =k + kiS + k 2 5 2 + 0(l)5 3 . 

3.2 Discrete Cole-Hopf transformation 

As observed by Lax [Tl], ^5], for the particular flux function / in (|3.1[) . 
one can perform a "discrete Cole-Hopf transformation which linearizes the 
scheme ()1.6j) . This enables us to write down explicit solutions of the Go- 
dunov scheme in the form of nonlinear superpositions of discrete traveling 
waves. More precisely, consider the Godunov scheme in this case, i.e. 



«r 1 = u i 



In [v + ^) - In [u + /ie u "-i) , (3.6) 



where we recall that [i > 1/2 and that v = 1 — fj,. We then have the following 
result due to Lax. 

Lemma 3.1. If z™ > is a solution of 



z 



n+1 - + vz*- , (3.7) 



j 



then a solution of (|3.6|) is provided by 
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Proof. A direct computation yields 



3 




z 



'J- 1 , i J z i-i 



= In + In 

V Z J J \v + fi^ 

= itf +ln (v + ^e"?- 1 ) - In [v + /xe n ^ 

□ 

In this section we will give discrete traveling waves for the the linearized 
scheme (|3.7j) and use these to construct discrete traveling waves for (|3,6j) . 
In the next section we will then use these as building blocks to write down an 
explicit solution of (|3.6j) which is a good approximation to the exact inviscid 
solution described in Section I3TT1 Letting both variables in 1)3. 7 Jl range over 
M we get 

z(t + 1, x) = fiz(t, x — 1) + uz(t, x). (3-8) 

Special solutions of (|3.8|) are easy to find. For every b > define a(b) as the 
Rankine-Hugoniot speed of the shock connecting the left state b to the right 
state u + = 0, i.e. 

a{b) := Wy+EH . (3.9) 
Notice that a = a(b) then satisfies the equation 

v + /J> = e^ b \ (3.10) 
and it follows that for any x\ G R the function 

Z(t,x) = e -bl*-xi-eiW 

is a solution of the linearized equation (]3.8[) . More generally, by linearity, 
any finite linear combination 

N 
3=0 

or any integral combination 

z(t,x) = [ e -c(*-*(0-<KC)t) dCj ( 3. n) 
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are also solutions of Q3.8|) . In turn, the formula 

"('■*'- ta ( £ W 2 ) (3 ' 12) 

allows us to recover an explicit solution of the nonlinear difference equation 



u(t + l, x) = u{t,x)- In (v + iie u{t > x) ) - In (v + fie uit ' x ~ 1] ) . (3.13) 

As an example let b\ < b 2 , si, X2 £ R and consider the solution 

:= e -6i[as-a;i-(T(6i)t] _|_ e -62[x-SB2-<T(6a)t] 

of l]3.8[) . A simple calculation shows that the corresponding function u in 
()3.12j) represents a traveling wave solution of (|3,13|) with speed 

, = ln(l + e bl )-ln(l + e fe2 ) = hafa) - b 2 a(b 2 ) 
h -b 2 h - b 2 

and connecting u{— 00) = b 2 to u{+co) =b\. To see where the center of the 
wave is located, we observe that at the point 

x x bi - x 2 b 2 61 

x = — = x 2 + — (xi - x 2 ) 

bi - 6 2 61 - b 2 

one has 

_ / e 6i e -6i(*-a!i) + e 62 e -62(«-«2)\ / e fei +e fe\ 6 1 + 6 2 

U(0, Xj = In r-rr r r— rz r = In 

v ' ' 1 g-fei(x-xi) _|_ e -b2(x-X2) I 

We can thus think of the point x*(t) := x + cr*t as the center of the discrete 
traveling profile at time t. 

4 Construction of a special discrete approximation 

We proceed to use the information from Section 13.21 to construct an exact 
solution of the upwind scheme (|3.6|) that is a good approximation of the 
particular solution of (|3,lj) with a single shock along the curve 7. This exact 
solution will be obtained as a solution u d (x,t) of the difference equation 
(|3.13j) as in ipnSj) . where z d (t, x) is a special solution of the linearized dif- 
ference equation (|3.8j) . (Superscript d is used (only) here to emphasize that 
the solutions satisfy the everywhere defined difference equations (|3.13|) and 
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(|3.8|1 .) We point out that the relevant property of the "good" approximation 
u d , which will be needed in the final computations, is that its level curves 
should travel with speeds sufficiently close to the level curves of the exact 
inviscid solution of Q3.1|) constructed in Section 13.11 

To construct such a solution of (|3.13|) we consider the curve j(t) given 
by (|3,2|) for large negative times. As above we set U-(t) := </?(7(i)), such 
that ct(u_(t)) = 7(7"), and define 

z(t,x):=l + J exp{-u_(r) [x-7(r)-7(r)(t-r)]}"0(r)dr, (4.1) 

for t G R. By (|3.11|) this defines a smooth solution of the difference equation 
1)3.81) . The corresponding solution of ()3.13j) will be denoted by u(t,x) (see 
Q4.14J1 below). For reasons that will be clear from the following computations 
we choose 

^(r) : = (-r)- 3 / 4 , r < 0. (4.2) 



To avoid several minus signs we let s = —t, 77 = —x, and set 



Z(s,rj) := z(t,x) 

[-2T 

1 + 



(4.3) 



exp 



[ u - (- T ) [v + l(- T ) ~ 7(- T )('S - T )] ]^(- T ) dr 
where we have made the change of variables r \—* — r. We define 
T(r) : 
W(t) : 

g(r) := ^(-r) = t~ 3 / 4 



- 7 (-r) = A r + 2^7, 
1 



u_(-r) = tp ( A + -j= ) , 



(4.4) 
(4.5) 
(4.6) 



for r > 0. Letting y := r] — T(s) and making the change of variables £ = r/s, 
we obtain 



Z(a,r(s) + y) 

r>2T 

1 + 



cxp 



[W(r) 



1 + 



W(r) 



y + T(s)-T(T)-t(T)(s-T 

y + v^(2- 



l + s 



1/4 



2T/s 



expi^(s£) [y + \/i(2- 



1 



I g{r) dr 
g(r) dr 

9(Qd£. (4.7) 
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We will consider times t in [— T, 0], i.e. s £ [0, T], such that the upper 
limit of integration in (|4.7|) is bounded away from £ = 1. The factor s 1 / 4 
will be important in the following computations and this is what dictates the 
particular choice in (|4.2I) . We proceed to analyze in detail the level curves 
of z(t, x) and u{t,x) by considering the level curves of Z(s,rj). For s» 1, 
the main contribution to the integral in (|4,7|) is provided within the region 
where £ ~ 1. A careful analysis is required in order to obtain the exact 
behavior. 



4.1 Level curves of u(t, x) 

As a first step we show that the time derivatives of the level curves of Z(s, rj) 
behave (almost) like T(s). Before starting to estimate the various terms that 
are needed, we make a simplifying observation. For e > 0, e <C 1, we let 
I(s) denote the interval 



I(s) = [l-s- Q ,l + s- Q ], 



a = 1/4 - e. 



Observation 1. In what follows all ^-integrals over [0, 2T/s\ \ I(s), as 
well as all boundary terms (obtained by differentiation of such integrals), 
are exponentially small. Namely, they are all of order 0(e~ CsE ) for some 
C,e>0. 



Next we define the functions 



Q{s,y) :-- 
H(s,0 :-- 



W(sO 



y + V^(2- 7e 



■ «i . 



Thus 



(e-i) 2 . 



f 2T/s 

Z(s, T(s) + y) = 1 + s 1/4 / exp s, y)] 5 (0 d£. 

J o 

Lemma 4.1. For y = O(l) and s>l we have 



(4.8) 
(4.9) 
(4.10) 



(4.11) 
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Proof. We start by stating a detailed expansion of the function F — Q — H 
where F, Q and H were defined by (jUJ), $M and <|^TU|) : 

F(e, S ,y)-Q( S ,y)- J ff( S ,0 = A + B(£ - 1) + C(£ - l) 2 + - l) 3 

+^-i) 4 + o(i)(e-i) 5 , 

where the coefficients A-i? are given as 

Mv,*) = y[— + 



s s 



3/2 



/ K l K 2 0(1) 

B(y,s) = -y[—L = + -± + 



lyfl S S 3/2 

. ki 3Kiy-2K 2 K 2 y + 0(1) 0(l)y 

n / \ _ Kpyfs 3«i 8k 2 - 5j/«i O(l) - K 2 j/ Q(l)y 
" 4 + 8 + 16Vi s + S 3/2 > 

Vs 

The proof is a direct calculation which is omitted. The estimate (|4.11|) is 
now obtained by exploiting the fact that, when differentiating with respect 
to time, we can subtract off the leading order contribution, i.e. the term 



d_ 

ds 



,1/4/ e ^?) d a 

Jl(s) I 



since this is exponentially small. We have 

Z(s, T(s) +y) = l + e ^V/4 f exp[F(e, s, y) - Q(s, y)]g(0 d£. 

Jl(s) 



T(s) 

such that, up to exponentially small terms, 

d 
ds 



Z(s,T(s)+y) = Q.eQ-sWf exp[F - Q]g(£) d* 

Jl(s) 

[ exp[F - Q]g(€) - e* <%] 

Jl(s) J 



ds 

o(i) , O(l) nl 



,3/2-e 



+ o(i) _ sl/4 /• (exp[F _ Q _^ b(e) _ 1)e //^ 

S J 1(a) 

+0(l)s 1 / 4 [ {{exp[F-Q-H}g(0-l)e H } d£, (4.12) 
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Using the expansion above and <?(£) = 1 — 3(£ — l)/4 + 0(1) (^ — l) 2 , we get 
that (for y = 0(1) and s » 1) 

expLF - Q - fT]5(0 - 1 = + {° dd function of (f - 1)}. (4.13) 

For the last term in (|4.12f) we apply the expansion above again to get that 
{(e W [F-Q-H]g(0-l)e H } s 

= [(exp[F -Q- H]g(0 -l)H s + exp[F -Q- H]g(0(F - Q - H) 

0(1) 



{odd function of (£ — 1)} + 



s 3/2-e 

Here we have used the expansion of g(£), 1)4. 13j) . that 

ko(^-I) 2 0(1) 
H s = — — = — — on I(s), 



and finally that 



(F-Q-H) s = ^-lf+ 



which follows from the expansion above. As 

«V4 



/ e ^)de = 0(l), 



r( s ) 

we therefore have that the two last terms in Q4.12|) are both 0(l)/s 3 / 2 ~ £ . 
The result follows. □ 

We now define 



A similar calculation as above shows that 

Z(s,r(s)+y) = l + c e y ^ + 



Q(s,y) , C^ 1 ) 



jl/2-e ' 

where 

/oo / 
exp^Ct, 0]« = 24/—. 
00 V K o 
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It follows that 



, ,'l + c e K °\ 0(1) Oil) 

U(s, T(s)) = In — ^— + =: a + 



1 + C J «V2-6 ' U s l/2- £ - 

For a sufficiently close to ao, and well within ] U(—oo), U(+oo) [, we denote 
the a-level curve of U(s,rj) by X(s;a), i.e. 

C/(s,X(s;a)) = a, (4.15) 

such that 

C/ S (s, X(s; a)) = -X s (s; a)U v {s, X{s; a)). (4.16) 
A computation shows that 

Z r ,( S ,T(s)+y) = K c e^ + ^^, 

whence 

K coe K °y(e KO - 1) Oil) 



U v {s,T(s)+y) 



(1 + K c e K oy) (1 + c e K os/) s i/2-e ' 

which is nonvanishing and O(l) for y = O(l). For a fixed s and with 
a = U(s, r(s) + y) it follows that 

±U(s,T(s)+y) = U v (s,T(s)+y)[X s (s;a)-t(s)] 

= 0(l)[X s (s;a)-t(s)], (4.17) 
for y = 0(1). On the other hand, using (|4.14jl and Lemma 14. 11 we get 
<*rr/_T^ , ^^( S ,r( S )+y + l) ^Z( S ,r( S )+y) _ O(l) 



-c/( s ,r( s ) + y) 



ds w y; Z(s,r(s) + y + l) Z(s,r(s) + y) s 3 / 2 ^ ' 

such that 

X s ( S ;a)-f( S ) = ^. (4.18) 

Transforming back to t and x coordinates, it follows that the a-level curve 
x(t; a) of u(t, x) satisfies 

xt(t; a) =7(t) + -ip_, (4.19) 

for a sufficiently close to ao- We conclude that there exists a smooth function 
c(o) defined on an interval about ao such that 

x(t; a) = c(o) + 7 (*) + (J^/L ' ( 42 °) 

for t G [-2T,0]. 
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5 Estimates on the total variation 



We now fix the solution u(t,x) := U(—t, —x) defined in ()4.14|) and set 

u] := u(n,j). 

By Lemma fH.ll these values provide a solution to the upwind scheme (|l.fi)) for 
the particular flux function f(u) = ln(z^ + fj,e u ). Recalling that we consider 
negative times we let v™ denote the corresponding solution of (|1.7jl obtained 
by prescribing vanishing u-data at time — 2T, and our goal is to show that 
the total variation of V(j) := i>j +1 over the set {j G Z | — T < j < 0} is of 
order lnT. Here T is as in Section^] This will be accomplished in a series of 
lemmas which reduce the problem to an estimate involving the heat kernel 
and the curve 7(t). In the course of doing this we will specify the precise 
assumptions on the function g(u) in (|1.7|) , 

Recalling that we have fixed A = 1/2 in (|1,7[) we get that the discrete 
Green kernel for (the linear part of) (|1,7|) is given by 

With vanishing data at time — 2T we get that the solution at time step n = 1 
is given as 

2T 

n=0 keZ 

where we have set 

/•u(n,k) 

r k ■■= -HO ~ 9{ul^)] =- / g'(£)dH. 

J u(n,k—l) 

We will make use of the following representation formula. 

Lemma 5.1. Let x(t;^) denote the ^-level curve of the particular solution 
u(t, x) specified above, and assume that the function g(u) in fl 1. 7| ) has com- 
pact support. We then have the representation 

poo r 

V(j)= / ^[E^i-I-C-nsOlJ*- 
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V(j) 



Proof. Using that the constructed solution u(t, x) is monotone decreasing 
at each time we get that, 

2T 




u{— n,k— 1) 



ru{— n,+oo) 

E 

27' 



j-[x(-n;£)] 



^)[E«?-W-n 



?1=0 



□ 



In the following computations we will use that the level curves x(t; £) of 
u(t,x) are well approximated by translates of 7(i), which is true for large 
(negative) values of i by l|4.2U|) . Now, for T sufficiently large, the values 
V(j) for j < —\[r are only influenced exponentially little by values of 
and x) for t G [— y/T, 0]. Furthermore, for j G [— T,— VT], 
the n-summations in 1)5.1 j) and 1)5.2)1 are unaffected, up to errors which are 
summable in j, by restricting n to the interval 



HI _ | 7 -|l/a+E \J1 + I -il/2+e" 



where < e <C 1 and /? = Ao — 1/2 > 0. This fact will be used repeatedly 
in the rest of this section. From now on all values j will be in [ — T, —y/T], 
i.e. far downstream. 

For two functions A(j) and B(j) defined for integers j we write A ~ B 
to indicate that the difference of the functions satisfy 

£ \A(j)-B(j)\<oo. 

j<-Vr 

With g as in Lemma 15. II we thus have 



V(j) 



V(o[ E 



j-[ir(-n;01 



ne/(j) 
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Since n G 7(jf) implies |2(j — [x(— n;£)J) — n\ < 0(l)n 1//2+e we can use 
Lemma 17. II in the appendix to conclude that 



V(j) - V(j 



f°° r 

« U „ T I A \ 



nel(j) 



l'(0[ E G x (%,2(j-lx(-n;0l)-n)}<% 

nel(j) 

Op -, 

</(£)[ E G xx (%,2(j-lx(-n;0l)-n)\<% 

nel(j) 
poo 

+ / g'(0 E o{n-^)di 

/(£)[ E G,(f,2(i-H-n;0])-n)]^ 

ne/(j) 

r°° r i 
" 4 / 9'(0[ E ^(|,2(j-[x(-n;0])-n)J^ 

ne/0) 
--:A(j)-B(j), (5.3) 



where we have used that = 0(\)\j\ l l 2+e and that n = 0{\j\) for n G 

We will estimate the variation of V(j) on [— T, — VT] by considering 
A(j) and £>(j) separately. 

5.1 Variation of B(j) 

We start by defining the functions 

nel(j) 
/•oo 

%;0 := / G xx (±,2(y-x(-t;0)-t)dt. 
Jo 

Lemma 5.2. VFe have 

H(j;0 ~ uniformly for £ in compacts. (5-6) 

Proof. We recall the notation /ii for the sawtooth function introduced in 



(5.4) 
(5.5) 
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ne/(j) 



Section [21 Using the decay properties of the heat kernel, we have 
\H(j;0-h(j;0\ = \ £ G xx {%,2(j-lx(-n;0l)-n) 

/ G M (|,20'-a;H;0)-t)t« 

JO 

£ [G w (§, 2(j - x(-n; £)) - n) + 0(l)G x . ra (§, n; £)) 
j) 

-/ G xx (±,2(j -x(-t;0)-t)dt 
Jo 

G xx {\, 2(j - x(-t; 0) - t) (1 + fci (*)) ^ 



nel(j) 



\G xxt {^2{j-x{-t-i))-t) 

+2[x t {-t;C) ~ \]G xxx {±,2(j - x (-t;0) - 



where n; £) is between 2(j — x(n; £)) — n and 2(j — [x(n; £)]) — n. Using 
the decay properties of the heat kernel, that xt = 0(1), and the fact that 
t = 0(l)\j\ for t G we conclude that 



\H(j;0 - h(j;0\ « 0(1)^ ^2+-j? dt 



Q(i) 

1^3/2- 



□ 



Note that the O(l) in the last term here is uniform with respect to £ as £ 
varies over a compact set. Since g is assumed to have compact support we 
conclude that the second term in (|5.3|) satisfies 



/(£)[ E G^(f,2(j-[*(n;£)])-n) 

ne/(j) 
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Using Q4.19JI we now calculate 

/>oo 

h(j;0 = I G xx (±,2(j -x(-t;£))-t)dt 

Gxx (§, 2(i - x(-t; 0) - t) [2x t (-t; - 1] tft 

/■0C 

+2 / G ra (|,20--xHO)-t)[l-i((-t;0]* 
•/ o 

poo 

/ C (|, 2(j - x(-t; 0) - t) [2xi(-t; - 1] 
•/ o 

+±G xt (t,2(j-x(-t;0)-t) dt 



OC 



+2 / ^(1,2(7 -x(-t;0)-t)((l-Ao)-^ + ^)(« 



Gxxx(|,2(j -x{-t;t))-t)dt, 



where we have added and subtracted G x t = G xxx . The first integral on the 
right vanishes so that rearranging gives, 

(2A -1)Mj;6 « "2 I™ G xx {±,2(j-x(-t;0)-t)^=+ [ ^ dt 

« d(= OW. (5.7) 



,m lil 2 UP' 2 -' 

Hence h(j;£) is summable in j, uniformly for £ in compacts. It follows 
from Lemma 15.21 that the same is true for H(j; £), and we conclude that the 
function B(J) is summable. Hence it contributes only a finite amount to the 
total variation of V(j) independently of T. 

5.2 Variation of A(j) 

It remains to show that the variation of A(j) in (|5.3j) grows indefinitely 
with T. This requires a more detailed analysis and we start by defining the 
functions 

Hvit) ■= ^^(§,2(7/ - x (-n;0)-n), (5.8) 
n>0 
/■oo 

:= / (?*(!, 2(y -a;(-t; (5-9) 

JO 

Note that the arguments in neither of these functions involve integer parts. 
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Lemma 5.3. We have 

~ l(j] Cj uniformly for £ in compacts. 



(5.10) 



Proof. As above we use the function h\ introduced in Section [21 and the 
decay properties of the heat kernel: 

/■oo 

\L(j; - Kr, 0\ = / G x (|,2(7 - x(-t; 0) ~t)(l + K(t)) alt 
Jo 

/■oo 

- G x (l2(j-x(-t;0)-t)dt 
Jo 

\G xt (^2(j-x{-t;t))-t) 

+ [2x t (-t; - l]G xx (|, 2(j - x(-t; £)) ~ t)] h^t) dt 

0(1) 



< 1 — ~ h 



/•oo 

/ [2x t (-t; £) - l]G xx (f , 2(j - 0) - t)M*) dt 
Jo 

To estimate the last integral we apply (|4.19f) and the function h 2 introduced 
in Section [21 to get 

/•oo 

/ [2x t (-t;0 - l)G xx (±,2(j-x(-t;0) ~ t)h x (t)dt 
Jo 

= ^2A -1 + ^) G«(f,2(j-x(-U))-^i(t)A 



(2A -1) / G xx (l2(j-x(-t;0)-t)ti 2 (t)dt + 



0(1) 

/•oo 

= 0(1) / G xxt (±,2(j - x(-t;0) - t)h 2 (t) dt 
Jo 

+0(1) r G xxx (^2(j- x (-t-i))-t)h 2 (t)dt+ 

Jo 

0(1) 

and the conclusion follows. 

Lemma 5.4. l(j;(,) ~ uniformly for ^ in compacts. 



b1 3/2- 



□ 
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Proof. Using (|4,19|) we have 

/>oo 

UO = / G x (l,2(j-x(-t;0)-t)dt 
Jo 

poo 

= / G x (l2(j-x(-t;0)-t)[2x t {-t;0-l]dt 
Jo 

/•oo 

+2 / G x {^2(j-x{-t;t))-t)[l-x t {-t;0]dt 
Jo 

G x (|, 2(j - x(-t; £)) - t) [2x t (-t; £) - 1] 

+|G t (|,2(3-a;(-t;0)-t)<ft 



+2 G x (| , 2(j - *(-t; 0) - t) ((1 - Ac) - ^ + ^) 



where we have added and subtracted Gt = G xx and used (|5.5|l . The first in- 
tegral on the right-hand side here vanishes. Recalling the definition of £), 
rearranging, and using the decay properties of the heat kernel together with 
what we have already proved about h(j;£), we get 



^fc.(Me-.(-«»-')|. 
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dt 



To estimate this last integral we repeat the same procedure: 

f°° rlt 
J G x [$,2(j-xi-t;t))-t)2= 

= j™ G x (1 , 2(i - x(-t; 0) - t) dt 
+2 G x (|, 2(j - z(-t; 0) - 

= {[2s t (-t;£) - 11^(1, 2(j - x(-t;0) ~ t) 

+ iG t (l2(j-x(-t;S))-t)} 

f°° rlt 

+2(1 -Ao) ^ ^(^(j-xH;?))-^ 

+0(1) / G s (|,2(i-x(-t;0)-t)? 
Jo 1 

1 f 00 rlt 

--] G xx {^2{j-x{-t-0)-t)^. 

Integration by parts and decay of G(t) shows that the first integral here 
is O ( 1 )/ 1 j 1 3 / 2 ~ e . Rearranging and using the decay properties of the heat 
kernel thus gives 

J r G .(i, a( ,-^ it ))- ( )*. j ^. 

We conclude that 

ittO- 0(1) 



|j|3/2-e ' 

which completes the proof. □ 

From Lemma f5.3l and Lemma l5.4l it follows that L(J; £) is summable with 
respect to j, uniformly for £ in compacts. We conclude that the function 

M(j;0-= E G x (%,2(j-lx(-n;0})-n) (5.11) 

nel(j) 
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appearing in the first term in (|5.3|) . may be expanded as 

n>0 
n>0 

where ((a)) := a — [a] denotes fractional part. Assume now that Ao is a 
rational number, 

\ _P 

*0 — ~> 

say, and write n = mq + i for 1 < i < q, to get that 

m>0 8=1 

•(x(-mg - 

Taylor expanding G xx about the points 

(ij 2 , 2 (j - x(-mq; £)) - mq) , 
and using the decay properties of the heat kernel, yields 



M(i;0 « ^ G **(^> 2 0" 0)-™?) 



m>0 



^E (^(-mg 



Recalling from 1)4.20(1 that the level curves x(t;£) are given as 
x(-t;0 = c{0-^ot + 2Vt + ^j^, fori>0, 



we obtain 

M(j;0 ~ £G^^2(j-x(-mg;£))-mg) 

m>0 

= ^G xx (^,2(j - x(-mq-,£)) -mq) 



i=i 



m>0 



53 (C(0 - f + 2y^ 



+ 1 + 



Q(ih 

Jmq), 



i=l 



(5.12) 
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Lemma 5.5. We have 



m>0 

uniformly for £ in compacts. 



E tto 



f + 2 \/™z) 



4=1 



Proof. By (j5.12[) it suffices to show that 
E | E 2 0' - *(-™*5 0) - mg) 

j m>0 

E «0 - 1 + 2 v^ + *K 0) - (c(0 - f + 2 Vm?) 

i=l 



< oo, 



where the summations are over those m and i for which mg + i 6 /(j). We 
have also set 

0(1) 



5(m, i) := 2\Jraq + i — + 



which is of magnitude 0(l)/yj7I f° r an m an d z in question. Next, for fixed 
£ and i, let 

and define 



a(s) := b + 2^ms, /3(s) := 6 + + £(s,i). 



Now, as s ranges over both functions and ((/?(s))) jump at most 

an 0(l)|j| e number of times. Also, if neither of these functions jump in an 
interval (s — l,s + 1), then they are 0(l)/y/\j\ close on that interval. It 
follows that the expression above is dominated by a sum of the form 



E 



\ o{i)W , Q(i)lil 1/2 



Y\ l J 'l 3/2 lil 3/ W 

which is finite. 

Making use of the identity 



□ 



i=l 



>-?] 



M + ¥ 
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and what we have already shown for H(j;£), we infer that 
M(j;0 ~ 52G xx (?X,2(j-x(-mq;Z))-mq) 



m>0 



{qc(0 + 2qy/mq) + 2=1 



^2G XX (^, 2(j - i(-m?; 0) " ™q) (qc(£) + 2 9v ^g) . 



m>0 



We next want to approximate M(j;£) with an integral and we proceed as 
above for H(j;£) and h(j;£). We define 

/•OO 

iV(i;0:= / G xx {f,2(j-x(-qs;0)-qs){qc(0 + 2q^)ds. 
Jo 

Lemma 5.6. M(j;£) ~ N(j;^) uniformly for £ in compacts. 
Proof. Using the function hi and decay properties of the heat kernel we have 
M(j;0-N(j;0 

/•OO 

« / G xx (3£,2(j-x(-qs;0)-qs){qc(0 + 2qVQt))h' 1 (s)ds 
Jo 

f°° d r i 
= - jf ^ [G M (f , 2(j - x(-qs; £)) ~ qs) (<?c(£) + 2q^)\h 1 (s) ds 

G xxx ( q i,2(j-x(-qs;0) ~ qs)[2x t (-qs;0 ~ 1] 
+ ±G xxt (f,2(j-x(-qs;0) ~ qs)} (?c(£) + 2q^ql) h^s) ds 
- jH G xx (f , 2(j - x(-qs; £)) - qs) { ^ (gc(£) + 2^)) } &i (*) & 

« - (f , 2(j - x(-g S ; 0) - ?*) { - £ rfs 

• y ° fc>0 

~ £ Gxxl^^O' -x(-gsfe;0) - QSk)hi{s k ), 
{k>o\s k ei(j)} 

where the s k denote the s- values for which qc(£) + 2q^fqs is an integer. 
Clearly s k = 0(k 2 ), such that the last sum contains only 0(l)|j| 1 / 4+£ terms. 
Since each term in this sum is 0(|j| -3 / 2 ) the conclusion follows. □ 
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Making a change of variables (r = qs) and restricting the integration to I(J) 
we get that 

M(i;0~-/ C(§,2(j-x(-T;0)-T)((gc(0 + 2(7v^))«!T. (5.13) 
Recalling that the level curve x(— r;£) is given as 

x{-r- = c(0 - A r + 2y^ + (5.14) 
we make a further change of variables r i— > cr where 



A straightforward calculation then yields 

G„(5,2(7 " x(-r;0) - r) = ^(2a 2 - l)e^ 2 + °W , 

which holds for r 6 ^(j), or, equivalently, for |<r — 2\/2| < 0(l)|j| £ . Here, 
and below, C denotes various explicit numerical constants. Substituting into 
()5.13)) we conclude that M(j;£) satisfies 



y^~1 /*00 

M(j;0 « ^ / (2a 2 -l)e" CT2 

1.71 J —CO 



*0 + Wf^ + l$S=J*' 

" W /-J 2 " 2 " 1)e " S (H + ^ + ^)) *■ 

We can now finally proceed to estimate the variation of V(j). We recall 
that the term B(j) in 1)5.3 J) is summable in j. Using the expression above 
for M(j; £) we can estimate the variation of the term A(j) in ()5.3j) . Setting 

Zj := 2g^/f (5.15) 

we thus get that 

/>oo 

u(j) - y(i - 1) « 4 / g '(OMU;OdC 

Jo 

= w\ ij 2(j2 - i)e ^ 2 ir 9 ' { ° i qc{ ° + Z] + "4 ^ 

C /*°° 2 

= — / (2cr 2 - l)e" a2 vr(cr; Zj ) da . 

\J I J — oo 
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Here 

z) := r g'(0 ((gc(0 + z + 3 ^ L } « , (5-16) 

*/ 

with c(£) as in 1)5.14)1 . Denning the function 

/oo 
(2a 2 - l)e- a2 Tr(a;z)da , (5.17) 
-oo 

we observe that U(z) is periodic with period 1 and (since Z-j = Zj) that 

Tot.Var. V(j) > V ^HM _ Cq , (5.18) 
-T<j<-VT J-^ j 

for some finite constant Cq independent of T. 

In order to complete the argument showing that V has large total vari- 
ation we need to introduce the following condition. 

Assumption (A): The smooth function g has compact support and is such 
that the 1-periodic function H(z) at {5. 17]) does not vanish identically. 

It is not difficult to see that this condition is satisfied for a large class of 
functions g. For example, if g' is formally replaced by a Dirac delta-function 
5 U * concentrated at a point u* close to ao, then we can directly compute 

H(z) = J °°j2a 2 - l)e~° 2 ((gcK) + z + fc^)) da £ 0. 

By continuity, any function g whose derivative approximates 5 U * will still 
satisfy the Assumption (A). 

We now fix a smooth function g which satisfies (A), and observe that 
the resulting function II is continuous. Therefore, for a sufficiently small 
number e, < e < sup 2 |II(z)|, there exist numbers a < b in (0, 1) such that 

|n(z)| > £ for all z G (n + a, n + b) =: J n , 

for all n > 0. Recalling the definition 1)5.1 5|) of Zj we define, for each n G N, 
the integers 

j(n) := smallest integer j for which Zj £ J n , 

and 

k(n) := largest integer j for which Zj G J n • 
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That is, 



3(n) 



~(3{n + a) 
4g 2 



We conclude from ()5.18|) that 

|n(zj)| 



E 



> 



Vr<j<r 



J 



> 



> 



and k{n) 



P{n + b)< 
4g 2 



|n(^; 



E E 

{n|J n C[^,T]}{il*j6J»} J 

e(k(n) - j(n)) 



E 

H J n c[^,T]} 

E 

{n| J„C[VT,T]} 



C 



n 



k(n) 



ClnT. 



This completes the proof and shows that there is no a priori bound available 
on the total variation of an approximate solution computed using the upwind 
scheme. 



6 Concluding remarks 

The analysis shows that it is possible to prescribe data for a strictly hyper- 
bolic system of conservation laws in such a way that the total variation of 
the solution generated by the Godunov scheme increases by an arbitrarily 
large amount, after a large number of time steps. We observe that the total 
change in the speed of the curve 7 decreases to zero as T — > 00. There- 
fore, as initial data for the u-component, we can take an arbitrarily small 
perturbation of a single shock. 

We summarize here the main features of our example: 

1. Our 2x2 system is strictly hyperbolic. One characteristic field is linearly 
degenerate, the other is genuinely nonlinear. The characteristic speeds are 
both strictly contained inside the interval [0, 1], so that the usual linearized 
stability conditions hold. 

2. We can choose a sequence of initial data of the form 

if x > , 
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where the functions <j) v are smooth and satisfy 

Tot. Var. ->• 0, ll<?V|lc fc -» 

for all k, as v —* oo. For these initial data, the exact solutions have uniformly 
bounded total variation, but the corresponding Godunov approximations 
(u u , v v ) satisfy 

Tot. Var. {y v (T v , •)} -> oo , 
for some sequence of times T;, — > oo. 

3. Consider a second sequence of Godunov solutions 

(u v , v v )(t,x) = (u u , v„)(t, x-1). 

By the previous estimates on the total variation, it trivially follows that 

\\v u (T u ) - v u (T u )\\ Tl 

hm —r. — —r. = oo . 

i'-oc ||u„(0) -u„(0)|| L1 

In other words, the Godunov approximations are also unstable in the L 1 
norm, w.r.t. perturbations of the initial data. This is in sharp contrast with 
the stability of Glimm or front tracking approximations ^E]) [Z|- 

Our counterexample has many special features: it can be solved in tri- 
angular form, the first equation contains a very particular flux function, 
and the second equation is linear. All these additional features allow us 
to perform explicit calculations and derive rigorous estimates, but none of 
them seems to be essential for the validity of the result. Indeed, we expect 
that the instability highlighted in this paper will be a common feature of all 
discrete schemes for systems of conservation laws, under generic conditions 
on the flux functions. In this direction, we observe that in our example the 
oscillations in the u-component are spread out on a very large spatial inter- 
val. Requiring that the first field be genuinely nonlinear (instead of linearly 
degenerate) would not achieve any appreciable decay in these oscillations. 

We observe also that our results do not show that the upwind scheme 
fails to converge. In light of the convergence results in 9 (obtained through 
compensated compactness arguments) we expect that one can in fact es- 
tablish strong convergence also for the systems considered above. In this 
connection it is worth noticing that triangular 2x2 systems of the form 
(|1.4|) - ()1.5[) are endowed with strictly convex entropies. 

A somewhat subtle point in the preceding analysis is the issue of small 
data. As presented above, our examples do not immediately fit the setting 



34 



of Glimm's theorem. We started by giving a particular shock-wave solution 
of the first equation (|1.4jl . and then chose a flux g(u) for the second equation. 
However it is now possible to go back and see that if we first fix a smooth 
function g(u), then we can give shock-wave solutions of (|1.4|) of arbitrarily 
small amplitude and with the properties: (i) its level curves travels with 
speed ~ [rational] + and (ii) the assumption (A) is satisfied. One can 

then redo the previous analysis for the resulting system/data pair. Since the 
data for the second equation are identically zero, this provides an example 
within the setting of Glimm's theorem for small BV data. 



7 Appendix 

Approximation of binomial coefficients For completeness we include 
the argument for the approximation of the discrete Green's function in terms 
of the heat kernel. Following the notation in Feller's book |lUj . 



a k (v) :-- 



2v 



2v 
u + k 



we have 



a k (v) 



a - 



exp 


E°=_ fc+1 rn(l-i); 


exp 


£? =1 i n (i-0; 





ao exp 



22i=l 3 



W 2 



Thus, 



a k {y) = a e » exp \ (£ ) - \ (f ) 



1 ffc^ 2 _ 1 fk\ 3 _ (fc-l) 2 fc 2 _|_ 1 (k}4 



4 \u 



Defining 



and recalling Stirling's formula 



tft(x) := 



-x 2 /2 



2ir 



139 



288n 2 51840n 3 



Oil) 



•nr 
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we get 



a k (u) = hm(kh)e £ 



Here e\ is due to the expansion above and £2 is due to the approximation 
in Stirling's formula. We have 



£1 = ei(k,u) :-- 



1 (k\ 2 1 fk\ 3 {k-lfk 2 „^,{k x4 



2\v ?>\v 



6v 3 



+ 0(1) 



while 



, . 1 1 O(l) 



Discrete Green's function Applying the approximation above to the 
discrete Green's function 



Kl := - 



(k — m)) e £l ~ £2 



we get for n = 2m, k = 21 that 
Kl = a k _ m (m) = 

= 2Gg,2fc-n) 
where G(t,x) = e~ x ' /2y/nt denotes the heat kernel and 



a £l-£2 



e 1 = £l (fc -!,!) = 



1 (2k-n\ 2 1 /2/c-n x3 



n 



n 



(2k-n-2) 2 (2k-nf 
12^3 



+ 



2k -n 



n 



and 



1 1 

£2=e 2 (f) = — - — + 



l2 4n 24n 3 n -1 
Lemma 7.1. 7/ \n - 2fc| < O (n l / 2+s ) (0 < 5 < 1), i/ien 

K£ - 7^ = 4[G X (§, 2fc - n) - G M (2 2fc - n) ] + O ( 



-2+4<P 
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Proof. Letting 

ei = ei(k -§,§), £2 = £ 2 (f), ei=ei(fc-n/2-l,n/2), 
we get that 



G(%,2k-n) -G(%,2k-n-2) 



3 £l-£2 



+2G (§, 2fc - n - 2) (e £l - £2 - e^ 2 ) 



= 4 



G x (§, 2fc-n) -G xx (§, 2fc-n) + 



O(l) 



£l-£2 



+2G(§,2A;-n-2) (e £1 -e fl )e _ea . 
Since (2/c — n)/n = O (n -1 / 2+<5 ) a calculation shows that 

ei - ^2 = O (n~^) , 

and that 



£l - 6! = O (n- 3 / 2+3S ) . 



Thus 



-^fc ~ #ra,fc-l 



G x {%,2k-n) -G xx (§,2fc-n) + 



O(l) 



(l + 0(n- 1+45 )) 



G x (2,2fc-n) -G xx (%,2k-n) + 0(1 



-2+45 



□ 
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